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The temperature and pressure dependence of the thermal displacements and lattice parameters 
were obtained across the phase transition of Ce using high-pressure, high-resolution neutron 

and synchrotron x-ray powder diffraction. The estimated vibrational entropy change per atom in 
the 7 ^ a phase transition, AS^^"' ~ (0.75±0.f 5)kB , is about half of the total entropy change. The 
bulk modulus follows a power-law pressure dependence which is well described using the framework 
of electron-phonon coupling. These results clearly demonstrate the importance of lattice vibrations, 
in addition to the spin and charge degrees of freedom, for a complete description of the 'y —> a phase 
transition in elemental Ce. 
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Materials with electrons near the boundary between 
itinerant and localized behavior continue to present a 
major theoretical challenge to a complete description of 
their properties, including multiple phases and anoma- 
lous thermodynamics. This is particularly true in the 4/ 
and 5/ systems, where this boundary appears to occur in 
or near the elements Co and Pu, respectively Q. In Pu, 
which possesses five allotropic phases at ambient pres- 
sure, a partial localization of some of the five 5/ electrons 
appears necessary to understand the higher temperature 
phases 0- Partial localization may also be present in 
U compounds j^]- Ce metal is in principle simpler, pos- 
sessing only a single 4/ electron, but still displays four 
different phases at ambient pressure. One of the most in- 
teresting and still not completely understood phenomena 
in Ce is the isostructural (fee) 7 — > a phase transition, 
which involves about 17% volume collapse at room tem- 
perature and pressure of roughly 0.8 GPa 0. 

In the majority of theoretical models 0, 0, IE S 
0, ^3 the 7 — > a transition has been attributed to an 
instability of the single 4/^ electron. The earliest mod- 
els focused on charge instability, while later models dealt 
with spin instability. The promotional model postulates 
a transition from Af^5d^6s^ (7-phasc) to 4/°5(i^6s^ (a- 
phase), but is inconsistent with the 4f binding energy 
and the cohesive energies of other 5(P6s^ materials [Sf. 
In the Mott transition (MT) model the 4/ elec- 
tron in the 7 phase is localized and non-binding, but is 
itinerant and binding in the lower volume a phase. The 
energy for the phase transition is provided by the kinetic 
energy of the itinerant / electron. In the Kondo- volume- 
collapse (KVC) model H 113 the 4/ electron is assumed 
to be localized in both the 7 and a phases, and the phase 
transition is driven by the Kondo spin fluctuation energy 



and entropy within the context of the single- impurity An- 
derson model. These early models ignored altogether 
an explicit treatment of the lattice degrees of freedom; 
even the lattice entropy is not considered. More recent 
treatments 0, include both the lattice and spin 
entropies, but still do not deal explicitly with the con- 
sequences of electron-lattice coupling despite the large 
volume collapse at the transition. 

In rare-earth compounds the electron-phonon cou- 
pling can be important because the ionic radii of dif- 
ferent valences often differ by over 10% For exam- 
ple, the anomalous phonon properties in mixed-valent 
Smo.75Yo.25S and YbInCu4 have been ascribed to a 
strong coupling to its valence instability 0, 0, 0, 0| . 
In addition, it has been shown that lattice vibrational 
contributions renormalize the two essential parameters 
of the Anderson model: the hybridization energy and 
energy of 4/ state 0, 0| . In 7-Ce, a comparison 
of phonon dispersion curves with those of thorium shows 
that the longitudinal branches of Ce are much softer than 
one would expect from the Lindemann rule |20j , which ac- 
counts for the differences in interatomic distance, mass, 
and melting temperature of these elements. This rel- 
ative softening, accompanied by a change in 4/ local- 
ization (see conclusions), suggests an important role for 
electron-phonon coupling, as in Smo.75Yo.25S, where the 
[CCC] longitudinal phonon branch is also soft compared 
to transverse branches [l^ . 

In this Letter, we further illuminate the role of the 
lattice in the 7 ^ a transition and present the first (re- 
markably) neutron diffraction experiment under pressure 
in elemental Ce. We measured thermal displacements of 
Ce in the 7 and a phases as a function of temperature at 
constant pressures (P ~ 0.4 and 0.5 GPa), and as a func- 
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tion of pressure at 300 K. The vibrational entropy change 
between the 7 and a phases, obtained from the thermal 
displacements using the Debye approximation, accounts 
for about half of the total entropy change and is thus 
non-trivial. Furthermore, the pressure dependences of 
the bulk modulus and thermal displacements are consis- 
tent with the explicit inclusion of electron-phonon cou- 
pling [Uli^l. These results suggest that lattice dynamics 
and electron-phonon coupling, in addition to electronic 
and spin instabilities, are essential ingredients for a com- 
plete understanding of the 7 ^ a phase transition in 
elemental Ce . 

A high-purity (99.99 %) polycrystalline, cylindrical in- 
got of Ce was obtained from the Ames Laboratory. In 
order to stabilize the 7 phase at room temperature, the 
Ce rod was sealed in a fused silica tube with low pres- 
sure (< 1 atm) Argon gas and placed in a furnace for 48 
hours at 423 K, and then furnace cooled over 6 hours. 
Powder diffraction patterns of the 7 and a phases were 
collected at room temperature in the pressure range up 
to 0.85 GPa on NPD at the Los Alamos Neutron Science 
Center (LANSCE) and up to 3.5 GPa on the X-7A beam- 
line at the National Synchrotron Light Source (NSLS), 
Brookhaven National Laboratory [23 • Temperature de- 
pendent measurements at constant pressure [P ^ 0.4 and 
0.5 GPa) were performed on SEPD at the Intense Pulsed 
Neutron Source (IPNS), Argonnc National Laboratory. 
For the neutron powder diffraction measurements, high 
pressure Al cells were used for the pressure range up to 
0.85 GPa. Contamination due to scattering from the 
cell was minimized using a Gd-based internal shielding, 
with spaces for the incoming, outgoing and ±90° scat- 
tered beam. Details of the cell designs are described 
elsewhere IMii]. 

Fig.^shows the isotropic thermal displacements (wf^o) 
of Ce as a function of pressure (upper panel) and tem- 
perature (lower panel) obtained from the Rietveld refine- 
ment using GSAS and EXPGUI The notable 
features in the upper panel are the very weak pressure 
dependence of (ufso) below 0.7 GPa and the pronounced 
drop of (wf^o) (« 30 — 40%) at the phase transition, 
which indicates a significant stiffening of the lattice. The 
square symbol is the thermal displacement calculated 
from measured phonon dispersion curves using the 
Born- von Karman force model at ambient conditions |28l| ; 
the agreement with our measurement is quite good. One 
finds that these measurements yield essentially the same 
relative change of thermal displacement (u^/u^) at the 
phase transition as the constant pressure measurements 
described next. In Fig. Q (lower panel) the closed squares 
represent measurements performed at 0.527 GPa with de- 
creasing temperature and open circles at 0.412 GPa with 
increasing temperature, respectively. At T=260 K (open 
circle) we had mixed 7 and a phases and were able to ob- 
tain the thermal displacements of both phases. The solid 
lines are fits using Eq.l^ [i^ with Debye temperatures 
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FIG. 1: (Upper panel) Isotropic thermal displacements of 
Ce vs pressure at 300 K. Symbols: circles are experimen- 
tal data; the square is calculated from Born-von Karman 
forces at ambient conditions figl . ^3 . Across the phase transi- 
tion the thermal displacement drops about 30%-40%. (Lower 
panel) Thermal displacements vs temperature at 0.527 GPa 
(closed squares) and 0.412 GPa (open circles). Data at 0.527 
GPa were collected with decreasing temperature and data at 
0.412 GPa were collected with increasing temperature. Data 
at T=260 K has a mixed 7 and a phases. Solid lines are 
fits to Eq. Q with Debye temperatures = 104(3) K and 
0£) = 133(3) K in 7 and q phases, respectively. 
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where $1 = j^"^'^ x{e^ — 1)^^ dx, Qd {—fu^o/kB) the 
Debye temperature, and (u^)offsct is a constant offset. 

At high temperatures, T > &d, the vibrational 
entropy can be approximated by ^vib ~ 3NkB[l + 
ln(T/8o) + ■■■], where 80 is the logarithmic phonon 
moment js^l- Approximating Qq hy Qd {Qd ~ 
Qqs^^^) yields the vibrational entropy change per atom, 
AS'Jj'j^" = ^vih ~ "^vib' which can be expressed as |3l| 



A C7-Q 
^"^vib 



3kB ines/e^. 



(2) 



01, 1^ / 133±3 K 



Using Eq. El we obtain AS^.^" 
(0.75 ± 0.15)fcB. This change is roughly half of the to- 
tal entropy change, AS'^ot" ~ 1.54A:b which follows 
from the latent heat or the Clausius-Clapeyron relation, 
dP/dT = AS'7„7"/AI/T-". This large change in vi- 
brational entropy is qualitatively consistent with sound 
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speed measurements [32, 123 in pure Ce (6-^ w 137 K 
and Oq « 154 K), but not with the measurements of 
the phonon density of states of the Ceg.gTho.i alloy at 
150 K (7-phase) and at 140 K (a-phase) by Manley et 
al. li^l, which showed little difference between 7 and a 
phases. We suspect that thorium atoms in Ceo.gTho.i 
alloys modify the elastic properties of Ce, as is observed 
in Mg-doped Ce alloys, where a few percent Mg signifi- 
cantly affect the Young's modulus, particularly in the a 
phase j^. 

We now focus on the softening of the isothermal bulk 
modulus Bt with increasing pressure in the 7 phase. 
The upper panel of Fig. 13 presents the measured P-V 
data. The closed circles in the lower panel of Fig. El 
were obtained by finite differences of the P-V data, 
B = -VdP/dV w ~VAP/AV. Here the errors are es- 
timated using the errors in applied pressures and lattice 
parameters. Lacking an equation of state that spans the 
7 — *■ a transition, we fit the 7 phase P{V) data to a 
cubic polynomial for a more accurate description of Bt, 
shown as solid line in the upper panel of Fig. |21 The 
Bt{P) derived from this fit is shown as solid line in the 

lower panel of Fig. |^ The dashed line in the a phase 

is from Olsen et al. [Sg- We note that our results are 
in good agreement with the adiabatic bulk modulus Bs 
obtained from ultrasound measurements by Voronov et 
al. lU [3^ . (The ratio between the isothermal and isen- 
tropic bulk modulus is Bt/Bs=Cp/Cv — 1 for Ce at 
room temperature.) Fig.|21shows that Bt decreases with 
increasing pressure below P < Pc, and that its smooth 
extrapolation in the 7 phase vanishes at Pc, indicating 
an elastic instability of the lattice. The discontinuity of 
Bt in the region P Ki Pc reflects the phase instability 
and the first-order nature of the transition. 

The softening of the bulk modulus with increasing 
pressure in the 7 phase is directly related to the soft- 
ening of the Cii elastic constant. In a cubic lattice the 
bulk modulus is B = [3Cii - 4C* -f- P]/3, where P is the 
applied hydrostatic pressure and C* = (Cn — Ci2)/2 is a 
shear modulus |38l. C * is known to be pressure insensi- 
tive in the 7 phase j2ll Is^ . Issj . Therefore, the softening of 
the bulk modulus in the 7 phase is a direct consequence 
of the softening of Cn. This result is consistent with 
the softening of the measured longitudinal sound speed 
with increasing pressure in 7-Ce [s^, 113 , and is related 
to the softening of longitudinal phonons |l9j at ambient 
pressure mentioned above. 

The pressure dependence of the thermal displacement 
and Bt of Ce can be described within the model of order- 
parameter-strain field coupling, e.g., XeQ^, where e is 
the strain, Q is a scalar order parameter, and A is the 
coupling constant. Bergman and Halperin showed that 
in this model a continuous phase transition is preempted 
by an elastic instability, causing a renormalization of the 
bulk modulus |2l| . (In rare-earth compounds, the order 
parameter Q has been defined as a relative change of the 
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FIG. 2: (Upper panel) P-V room temperature isotherm. The 
solid line is a cubic polynomial fit to the 7 phase of Ce. (Lower 
panel) Bulk modulus vs pressure. The symbols are calculated 
from finite differences. The solid line is obtained from the 
polynomial fit to the 7 phase P-V data in upper panel. The 
dashed line in the a phase is from Olsen et al. |3fl |. Inset: 
The dashed line is a fit to the solid line near the transition 
(0.7 GPa< P < 0.837 GPa), where ^ \P - P^\^°' with 
Pc « 0.83 GPa and a ^ 0.46. 



lattice constant (Aa/a) [isl H^ . which is proportional 
to the 4/ occupation number n / [3^ ^3 ■ Because the 
strain couples to the square of the order parameter, it 
constitutes a "secondary" order parameter [131 ■) In this 
framework the thermal displacement behaves as (u^) ^ 
Ei+ E2p-\- Esp^'" , wherep =| P — Pc |, are constants 
and a is the specific heat critical exponent • The bulk 
modulus vanishes as B ^ p°'. In mean- field theory the 
exponent is a = 0, while for Gaussian fluctuations one 
expects a = \ [2^ . 

As shown in Figs. Q] and |21 our experimental data 
exhibit the features of such a model. Our pressure- 
dependent measurements at 300 K show no divergence of 
the thermal displacement, instead, they exhibit a sharp 
drop across the first-order transition. In addition, we 
find a continuous decrease of the bulk modulus towards 
zero as the pressure approaches Pc with a simple scaling 
behavior, Bt ~ (see inset in lower panel of Fig. [SJ. 
Near the transition we extract the exponent a w 0.46 
and the critical pressure Pc ~ 0.83 GPa. Of course, a is 
not a truly measured "critical exponent" because (1) we 
have only a few measured data points below the transi- 
tion, and (2) the descent of Bt towards zero is clearly 
preempted by the first-order phase transition, as seen by 
the large volume collapse. 
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In conclusion, these new experiments and analysis 
clearly demonstrate that vibrational entropy plays 
a significant role in stabilizing 7-Ce, accounting for 
about half of the total change in entropy at the 7 — > a 
transition. The pressure dependence of the thermal 
displacement and bulk modulus of Ce also strongly sug- 
gest that electron-phonon coupling plays an important 
role in the 7 — > a transition. Finally, we note that our 
results mean that a complete understanding of Ce and 
its fascinating 7 ^ a transition must take into account 
the important interactions between the spin, charge and 
lattice degrees of freedom, and that this relates closely 
to the competition between localization and itinerancy 
in /-electron materials. 
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